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Abstract 

The production of d' dibaryons in heavy ion collisions due to the elementary 
process NN d'lr is considered. The NN —> d'n cross section is estimated 
using the vacuum d' width T^/ ~ 0.5 MeV extracted from data on the double 
charge exchange reactions on nuclei. The d' production rate per single collision 
of heavy ions is estimated at an incident beam energy of 1 A-GeV within the 
framework of the Quantum Molecular Dynamics transport model. We suggest 
to analyse the invariant mass spectrum of the NNtt system in order to search 
for an abundance of events with the invariant mass of the d' dibaryon. The 
d' peak is found to exceed the statistical fluctuations of the background at a 
6(7 level for 2 10® • A central collisions of heavy ions with the atomic number 
A. 
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I. INTRODUCTION 


In the last years the possibility of the existence and the properties of a narrow dibaryon 
d' (T = 0, = 0“) have been widely discussed both from the experimental and theoretical 

points of views [0] - [§]• So far, the narrow resonance structure is unambiguously established 
only in the double charge exchange (DCE) reactions on nuclei from to |^. The 
peaks in the forward DCE reaction cross sections at the incident kinetic energy of pions 
Ttt ~ 50 MeV are well described assuming that the dominant contribution of this process 
stems from a d' dibaryon with a mass M^’ = 2063 MeV and a free decay width ^ 0.5 
MeV. The d' mechanism is presently the only one which allows to explain the observed 
narrow structures. It is therefore important to hnd an evidence for a d' dibaryon in other 
reactions. At CELSIUS (Uppsala), the reaction pp —>• ppn^TT~ was studied |@]. The analysis 
of the ppvT” invariant mass spectrum revealed a 4cr peak ai M k, M^'. The reported data are, 
however, still preliminary. Further plans to search for a d' dibaryon include proton-proton 
collisions and photo- and electro-production on deuterons. 

The existence of narrow dibaryon resonances has important consequences for physics of 
dense nuclear matter. When the density is increased beyond a critical value, the formation 
of a Bose condensate of the dibaryons becomes energetically favourable. This phenomenon 

- 0 - 


has been studied in Refs. [[ 

In the present work, the possibility of the d' production in heavy ion collisions is investi¬ 
gated for incident beam energies around 1 A-GeV. Such a beam is, e.g., available at SIS at 
GSI (Darmstadt). We investigate the d' production rate and the pp'K~ background within 
the framework of the Quantum Molecular Dynamics (QMD) transport approach |]^ which 


is well established and describes successfully many observables of heavy ion reactions 

g. 


In the next Sect., the W —>■ d'n cross section is calculated using as an input value the 
vacuum d' width Td' ^ 0.5 MeV extracted from the DGE reactions. We expect 2 10“'^ d' per 
single NN collision. In Sect. 3 we start with a qualitative discussion of the d' production 
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dynamics in heavy ion collisions. Using reasonable assnmptions, we give an estimate for the 
d! prodnction rate and for the nnmber of heavy ion collisions needed to separate the d! peak 
from the backgronnd. 

In Sect. 4, the d' prodnction in heavy ion collisions is stndied with the nse of the QMD 
simnlations. The d' prodnction rate is calcnlated for central collisions of the ions 

and Au. The invariant mass distribntion of the backgronnd pp'K~ events 
is analysed in Sect. 5. We found that it is possible to distinguish a 6 cr d' peak from the 
analysis of at least 2 10^ • A collisions. Corresponding experiments at the SIS have collected 
up to now about 10® events. Thus the possibility for the search for the d' dibaryon in heavy 
ion collisions is apparently a realistic aim. 

II. CROSS SECTION FOR D' PRODUCTION IN NUCLEON-NUCLEON 

COLLISIONS 

The possible NNd'Ti couplings are carefully analysed in Ref. 0. Thus we start with the 
form of the NNd'n vertex suggested in this work 

TNNd'n = gillC'^^Tii)2'Kid'. ( 1 ) 

Here, C = 7270 and 5 ^ is a form factor which depends on the momenta of the particles 
involved. The two nucleons in the vertex (0 are in the relative ^S'o state, so they are 
affected by the long-distance attraction resulting in a virtual pole and by the short-distance 
repulsion. The hnal-state interaction of two nucleons substantially modihes the d' decay rate. 
On the other hand, the initial-state interaction of the nucleons in the reaction NN ^ d'n 
modihes the cross section. 

First we assume a constant form factor g and neglect the interaction of the nucleons. 
Then we discuss in details the hnal- and initial-state interactions and the inhuence of the 
nontrivial momentum dependence of the form factor g. 

The total width T^?^ of the d' dibaryon for a constant g with no hnal-state interactions 
can be calculated as 
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( 2 ) 


ri°’ = 


QAir'^ 

Here, Tj^^max = {{Md> — — Am?) /(2Mrf/) 45 MeV is the maximum pion kinetic energy in 
the d! decay, p, m, and Md' are the pion, nucleon, and dibaryon masses, respectively. Due 
to the isotopic invariance, the partial decay widths are given by 


p(o) _ p(o) 

d'—>pp7r“ d'—^pmr' 


_p(o) _ 

0 ■*' d'—>-nmr+ 3^ ' 


The same relations hold true for the widths corrected to the hnal-state interactions. 
The total cross section of the dibaryon production equals 


a. 


(0) 


1 k 


pp^d'n+ 9 


( 3 ) 


where k and p are the pion and nucleon momenta in the center-of-mass frame. In the 
different isotopic channels, the d! production cross sections are related by 


r(0) 


_ ^(0) 


^pp^d'TT+ ^nn—fd'TT~ ‘^^pn^d'-K^' 


The same relations are valid for the cross sections corrected by the initial-state interactions. 
In terms of the d' width, the pp cross section takes the form 


/o) 

pp^d'iT^ 


yo) Svr 

d' 


k 


( 4 ) 


3 y/^^Md'Tl^^^p 

If the hnal-state interaction of the nucleons in the d' decay and the initial-state interaction 
of the nucleons in the d' production are taken into account, we obtain both, a corrected width 
Trf/ and a corrected cross section app^d'-K+'- 


Td' = and (Jpp^d'^+ = (^pJ^d'n+Vi- 


( 5 ) 


The corrected cross section contains, as compared to Eq. (^, an additional factor pi/pp 

Stt 1 1 k Pi 


>d'7r+ h, 


d'- 


( 6 ) 


3 v'Wf P')F 

In contrast to Ref. |0] where the factors pi and pp have been determined form different 
sources, here we use a method which allows to calculate both, pi and pp, consistently. The 
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common scale of rjj and rjp remains in our method unfixed due to the unknown asymptotic 
behaviour of the nucleon scattering phase shifts at high energies. The common scale factor, 
however, drops out from the ratio rji/r^p entering into Eq. (^). The present method is based 
on the integral representation of the dost function ||T7[ 


/ 1 roo 

Di{T) = exp- 


d(T') 


TT do T' — T 


dr . 


(7) 


Here, T is the nucleon kinetic energy in the laboratory frame. The value of T is connected 
to the invariant energy s = (pi + ^ 2 )^ by the relation s = AnT + 2mT. In Eq. (^, 5{T) is 
the phase shift in the dominant partial wave state of the two nucleons in the d' decay. 

It is known from scattering theory that the influence of the hnal-state interaction on the 
outgoing particles can be taken into account by dividing the bare amplitude by the dost 
function. The same prescription is valid for the initial-state interaction. The amplitudes 
with two external nucleon lines get therefore factors 1/Di{T). 

The enhancement factor rjj entering Eq. (|^) is simply 

“ |L)/(T)|2’ 

while the factor rjp can be calculated by averaging the same function ri{T) over the phase 
space distribution of two nucleons originating from the decay d' —>■ NNtt. 

The phase d(T) is experimentally known up to T = 1000 MeV |]T^ and can be extrap¬ 
olated to inhnity in a model dependent way. We use the experimental data in the interval 


T = O-h 1000 MeV and results of the model 19 for T > 1000 MeV. The uncertainties in the 


asymptotic behaviour of the phase do, however, not signihcantly influence the ratio r]i/rjp. 
The function (5{T) = \\n.{r]iT )is shown in Fig. 1 for the interval T = 0 -h 1000 
MeV (solid curve). It is seen that it rapidly falls off from /d(0) = 0 to a minimum value 
of fd{Ti) ^ —3.3 at a kinetic energy of Ti ^ 600 MeV. In the interval T = 700 -=- 1000 
MeV, (d{T) is a slowly varying function. For comparison, we plot in Fig. 1 the function 
/3(T) calculated in the effective-range approximation (see Appendix). The ratio rji/rjp in 
our approach becomes 
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( 8 ) 


Vi 1 
Vf 5' 

Inserting Md' = 2063 MeV and = 0.5 MeV, as suggested by the interpretation of the 
DCE reactions into Eq. (^) we obtain for the cross section 


Upp^d'nF 


k 

300— fib. 
P 


( 9 ) 


In the reaction NN —> d'vr, the threshold energy is Tq = 710 MeV. At T 730 MeV, 
and k/p k. 1/10 one obtains app^d'-K+ = 30 fih. The cross section (P) should be compared to 
the experimental data on the d' production in the pp pp7r~^F~ reaction P]. At T = 750 
MeV, where the total cross section of the pp —>• ppvr+Tr” reaction is about 10 -i- 20 fih, the 
d' contribution is estimated to be at a 7% level. It means that at this energy app^d'TT+ ~ 
15 ■ 0.07 -3^3 fih, i.e. the empirical value is about ten times smaller than the value given 
by Eq. (j^). Hence, we are forced to conclude that the constant form factor g does not 
allow to calculate the cross section app^d'n+ correctly from the width F^/. Thus, we need to 
introduce a form factor to suppress the NNd'n vertex at high NN energies. One possible 
choice of the form factor is the following 
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T2 

—2— 


( 10 ) 


where rriFF is a cut-off mass. The required suppression allows to determine the mass pa¬ 
rameter rriFF ^ 240 MeV. Finally, we get for the cross section 


>d'7r+ 


300 


(i + Ai) 


-If,*. 

'P 


( 11 ) 


Fig. 2 shows the cross section versus the nucleon kinetic energy from the d'lr threshold up 
to T = 1000 MeV. 

It is seen from Eq. (7) that the enhancement factor rf{T) approaches unity with increasing 
T. In the potential framework, the hnal- and initial-state interaction effects disappear with 
increasing the energy. In our case Vi/vf ~ 1/5 and rfi is almost constant in the interval T = 
700 -i-1000 MeV. One can thus assume that at T = 700 -i-1000 MeV initial-state interaction 
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effects are essentially switched off. The result rji/riF k. 1/5 should be compared to the result 
of Ref. 0: rjp 5, 77 / ~ 1/10, and rji/r]F ^ 1/50. The order-of-magnitude difference is 
apparently related to the effect of a d! form factor implicitly included in calculations of the 
enhancement factor rjj performed in Ref. |p. The effective-range approximation gives rjF ^ 3 
and rjj ^ 1 (see Appendix). 

III. NAIVE ESTIMATE FOR D' PRODUCTION IN HEAVY ION COLLISIONS 

In the previous section the cross section pp —dV"*" has been estimated using thereby 
as an input value the vacuum d' width T^/ 0.5 MeV. This cross section turned out to be 

rather small. For a kinetic energy of T = 1 GeV of the hrst nucleon in the rest frame of the 
second nucleon the d' production cross section reaches a value of app^d'n+ ~ 10 /ib (see Fig. 
2). The total nucleon-nucleon cross section at the same energy is about ~ 40 mb. The 
ratio between the total d' production cross section averaged over the isospin states of the 
initial nucleons and the total nucleon-nucleon cross section is about 


We thus expect ^ ^ 2 10“^ d' dibaryons to be produced per single NN collision. 

In heavy ion collisions the number of the elementary NN scatterings is already at in¬ 
termediate energies relatively large and thus the d' production rate is enhanced. The mean 
free path 1 fm of a nucleon in a heavy nucleus is small as compared to the nuclear size. 
Therefore a nucleon from the colliding nucleus undergoes several scatterings with nucleons 
from the target nucleus, loses its energy and hnally cannot produce a d' dibaryon (y^ < M^/) 
anymore. A naive estimate of the number of the d' dibaryons per a single AA collision yields 



(13) 


The d' production rate in symmetric and central heavy ion collisions is roughly enhanced by 
a factor A corresponding to the number of the projectile/tar get nucleons as compared to pp 
collisions. The background in heavy ion collisions is more complicated. 
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We wish to select NNn events to test the invariant mass distribution with respect to 
the presence of a peak caused by the d' dibaryon. For total number of the pions produced 
per single AA collision at an energy of 1 A-GeV, we assume N.,^ ~ A/d. This number is an 
upper limit for the pion production rate in the central heavy ion collisions, which coincides 
roughly with the prediction of the so-called Harris systematics at that energy. However, 
recent experiments indicate that pion production is suppressed in heavy systems with 
respect to the Harris systematics. The total number Nt of those NNtt events which have to 
be analysed is then 


/ A 

Nt = 2A^N^Nc ^ i 10^ i - Nc, 

V2OO/ ’ 


(14) 


with Nc being the number of nucleus-nucleus collisions. 

For the further discussion we denote with ^.^(p) and nB(k) the momentum space distri¬ 
bution functions of the outgoing nucleons and pions, respectively. The number density of 
NNn triplets with respect to the NNn invariant mass Mt is given by 


1 dNt 


9nucl9pion 

2 ! 


^^ ""(p') ""(p") 


(16) 


where gnuci = 2x2, Qpion = 3 are statistical factors for nucleons and pions. Since we are 
interested in the values of Mt close to the NNtt threshold, one can write 


1 dNt 


9nucl9pion 


2! (27r)' 




dP 


Mt 


(27r)3 + P2 


n-A^P) riBi-I^P) 




'Mt 


(16) 


where 


3 

$3(M,) = I dptdp2dkS\pt+p2 + k-P)^ ^{pMtf^Tl^,^{Mt) (17) 

is essentially the phase space of the NNtt system. Here, P = {\J Mt -|- P^, P), m and p are 
nucleon and pion masses, and T^r^max = — pY ~ 4m^)/(2Mi) is the maximum kinetic 

energy of the pions in the center-of-mass frame of the NNtt system for a given invariant 
mass Mt- 
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Near the NNtt threshold, the dependence of the density dNt/dMt on the invariant mass 
Mt is mainly governed by the factor One can therefore write 


1 dN, 




4 10““ MeV 


-1 


N, dM, Tl^(M,)dM, 

The density of the background events in vicinity of the d' peak is estimated as 

1 dNt _ 2A^N^ dNt _ 1600 / 
l^jMt ^ Nt dMt MeV ‘ 


(18) 


(19) 


In the case of the Au + Au collisions, we thus obtain about 1600 background events 
per 1 MeV bin of spectrum of the invariant mass Mf. 

The actual bin width AMt which we are interested in is given by the maximum of the 
d' width Trfi 0.5 MeV and the experimental resolution The number of the 

background events in this interval. 


ANt 2A^N^ dNt ^ f A f AMt \ 

- - = - -—-AMt ^ 1600 i -- 

Nt dMt V200; UMeV; 


N< 


c 


( 20 ) 


should be compared to the number of the d' events given by Eq. (13). 

Suppose we have Nc collisions and we want to have a statistically signihcant peak at a 
no-cr level with > 6. We should require 


fN NcAni < NcHd', 


( 21 ) 


with Arit = ANt/Nc being the number of the NNtt events per single AA collision in the 
interval AM*. We get 


Nc > n; 


Ant 


4 10^ 




/2 10-nV AMt \ 

\ ^ ) UMeVW 


( 22 ) 


The above estimates give a rough idea on the number of collisions needed to reveal the 
d' peak. We see that light nuclei require smaller number of the collisions and therefore are 
better candidates for d' searches. In order to obtain a more precise estimate and to check 
the A dependence of the results, we perform in Sects. 4 and 5 Monte Carlo simulations for 
the d' production in central heavy ion collisions within the QMD model. We analyse also 
additional cuts and show that they allow to reduce the value Nc by a factor 1/2. 


9 



IV. SIMULATION OF D' PRODUCTION IN HEAVY ION COLLISIONS 


The simulations were done within the framework of the Quantum Molecular Dynamics 
(QMD) model |^. QMD is a semi classical transport model which accounts for relevant 
quantum aspects like the Fermi motion of the nucleons, stochastic scattering processes in¬ 
cluding Pauli blocking in the hnal states, the creation and reabsorption of resonances and 
the particle production. The time evolution of the individual nucleons is thereby governed 
by the classical equations of motion 

dH dQi dH 


dt 


_^ ^_ 

9qi ’ dt 


(23) 


with the classical iV—body Hamiltonian H 

» = E VpTFw? + i E + U't + Upr') ■ 


(24) 




The Hamiltonian, Eq. (|^) , contains two-body potential interactions which are hnally de¬ 
termined as classical expectation values from the local Skyrme forces f/^ supplemented by 
a phenomenological momentum dependence. 


Uij = a 


Pij 

.Po, 


\P0, 


+ 5ln^ (e|pi - Pjp 


pij 
PO '' 


(25) 


and an effective Coulomb interaction 


jjcoul _ / ^ 

^ij - \ ^ 


e2 


Iq* - 


-erf 


|q* - cb'l 




(26) 


Here, pij is a two-body interaction density defined as 

1 


Pij 


(dvrL) 


,-(qi-qi) /4i 


(27) 


and erf is the error function. The parameters a,f3,'y,6,e of the Skyrme interaction, Eq. 
(p5|), are determined such as to reproduce the saturation density (po = 0.17 fm“^), the 
binding energy {Eb = —16 MeV) for normal nuclear matter at a hxed incompressibility 
Koo, as well as the correct momentum dependence of the real part of the nucleon-nucleus 
optical potential . In the present calculations we use a soft equation of state with = 
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200 MeV corresponding to the momentum-dependent interaction (SMD), Eq. (|25|) , which is 
known to reproduce reasonably well the reaction dynamics in the considered energy range. 
The Yukawa-type potential in Eq. (p^) serves to improve the surface properties and 
stability of the initialised nuclei. 

Two-body collisions are implemented by Monte Carlo methods, the Pauli principle is 
thereby taken into account for the hnal states. Nucleons while propagating collide with each 
other stochastically provided that the distance between centroids of the two Gaussian wave 
packets is less than dmin = \J'ytoti'/s)/'n. For the inelastic nucleon-nucleon channels, we 
include the relevant baryonic resonances, i.e. A(1232) and iV*(1440). The pions originate 
from the respective resonance decays. In particular we include the resonance production and 
rescatterings in the inelastic NN collisions: the one-pion decays of A and A*, the two-pion 
decays of A*, and the one-pion reabsorption processes. For details see Refs. [ 0 . 0 - 

Since the production cross section of the d' dibaryons is small, these particles are treated 
perturbatively. The d' dibaryons are produced at kinematically favourable conditions in the 
pp and nn collisions with a probability 


a 



P 


.tot 

NN 


while in the pn collisions the d' dibaryons are produced with a probability P/2. Here, a] 


is the total (elastic -I- inelastic) NN cross section and app^d'-K+ is given by Eq. (|TI]). The 
results of our simulations for the d' production rate in 1000, 50, 25, and 100 central collisions 
of the ions ^°Ae, ®®Ai, and Au are shown on Fig. 3 (a). The impact parameter b is 

assumed to be less, respectively, than 0.7, 1.3, 2, and 3 fm (minimal bias). In these intervals. 


the probability for the d' production is not sensitive to the values of b. The d' production 
rate increases linearly with the atomic number A, and so our results can be presented in the 


form 


Ud' 10 ■ A. 


(28) 
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V. THE INVARIANT MASS DISTRIBUTION OF THE BACKGROUND PPtt- 

EVENTS 


The number of the d' dibaryons produced per single collision is much smaller than the 
number of the NNtt background events in the vicinity of the d' mass (according to Sect. 3). 
It is therefore reasonable to search for cuts which are able to enhance the d' signal. 

The hrst possible cut takes the fact into account that the d' is an isosinglet state. Only 
1/12 of all NNti events form an isosinglet state (we assume isotropy of the NNir events in 
the isotopic space). The selection of isosinglet NNti events could reduce the background by 
a factor of 12. Such a selection is, however, impossible. Let us select only pp7i~ events. The 
total reduction of the background is also 1/12. Since the branching ratio of the d' pp7r~ 
decay is 1/3, the number of the d' events is reduced as well. The relative reduction of the 
number Nc of collisions needed to reveal the d' peak is, according to Eq. (|^), 3^/12 = 3/4. 

The second possible cut can be obtained by the selection of protons with an invariant 
mass Mpp < 2m + 15 MeV. Such a cut not much influences the d' signal, since for the 
dominant part of the d' decays (~ 2/3) the invariant mass of the two protons Mpp is less 
than 2m + 15 MeV. This is just the effect of the final-state interaction of two protons in the 
d' —>■ ppTi~ decay. In Fig. 1, this effect is displayed in the fast increase of the function /9(T) 
at small nucleon kinetic energies T. The reduction of the background can be estimated as 
in Eq. ([l6l) and gives 


T^{Mpp = 2m) - T^{Mpp = 2m + 15 MeV) _ 1 

T^{Mpp = 2m) ^3 ^ ^ 

where T^^^Mpp) is the kinetic energy of pions in the NNtt system for a given invariant mass 
Mpp of the two nucleons. In particular, T^^^Mpp = 2m) = T^r^max ~ 45 MeV. The relative 
reduction of the value Nc is (1/3) ■ (3/2)^ = 3/4. 

The d' dibaryon decay products are distributed isotropically in the nucleus-nucleus 
center-of-mass frame, while the nucleons in heavy ion collisions are produced with a higher 
probability in the forward and backward directions. The angular distribution of the outgoing 
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protons averaged over 100 Au + Au collisions is shown in Fig. 4. The cnt in the polar 
angle of the beam can provide an additional redaction of the lower limit for the nnmber of 
the collisions repaired to reveal the d' peak. The backgronnd is rednced by a factor 2, while 
the nnmber of the d' events is rednced by a factor 3/4 (see below). The relative redaction 
of the nnmber Nc is small (1/2) ■ (4/3)^ = 8/9. 

These three cuts when applied separately are not much efficient. Nevertheless, being 
combined they reduce the value Nc by a factor (3/4) ■ (3/4) • (8/9) = 1/2. 

The simulations of the ppTT~ background are made within the framework of the QMD 
model for 1000, 50, 25, and 100 central collisions of the ions ^^^Ag, and Au 

at the beam energy 1000 MeV per nucleon. In order to exclude complex fragments from 
the analyses, we take only protons into account which have no neighbouring nucleons closer 
than 3 fm in the hnal stage of the reaction. Here we trace the evolution of the reaction up to 
80 fm/c, a time where almost all resonances (A and A*(1440)) have decayed into nucleons 
and pions. 

The invariant mass distribution of the pp'K~ events in the interval 2.02 < Mt < 3 GeV 
is shown in Fig. 5 (a) for one typical Au collision. This distribution has a 

maximum at Mt ^ 2.3 GeV. The upper limit in the integral entering the denominator of 
Eq. (|T^ is therefore reasonable. 

In Fig. 5 (b) we show the invariant mass distribution of the events in the interval 
2.02 - 2.1 GeV, averaged over 100 collisions Au + Au. According to Eq. (^), the 
density of the NNtt events is proportional to the phase space of the NNtt system. Near the 
NNtt threshold, the phase space is proportional to square of the difference Mt — Mth where 
Mth = 2m + /i. We see that in this region the density of the ppTi~ events is htted well by the 
curve C{Mt — Mth)’^, in agreement with the expression ( 0 )- From the £t, we determine the 
coefficient C = 2.6 10^ GeV“^. In vicinity of the d' dibaryon mass, the average number of 
the background pp'K~ events equals ANt/Nc = 55. This number should be compared with 
the number in the right hand side of Eq. (j^). Here we got 1600 NNtt events per 1 MeV 
or, equivalently, 1600/12 = 130 pp7r~ events per 1 MeV. The estimate of Sect. 3 is therefore 
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in reasonable agreement with the results of the QMD simulations. 

In Fig. 5 (c) we show the distribution of the ppn~ events for 100 collisions Au + Au, 
grouped in 10 bins of a width of AM* = 1 MeV around the d' dibaryon mass. The solid 
curve is the best linear £t to this distribution. This distribution is now further analysed in 

Fig. 6. 

First, we construct the histogram (a): The number of the bins versus the deviation of 
the number of the ppTT~ events from the mean values determined by the solid straight line in 
Fig. 5 (c). In the vicinity of the d' mass, the number of NNtt triplets is about ANt ~ 5500. 
From Fig. 6 (a) we deduce the RMS deviation a ~ 80 in the number ANt. The Poisson 
law requires a ~ \/ANt ~ 75. The background fluctuations are therefore distributed in 
agreement with the Poisson statistics. 

In Fig. 6 (b), the same results are shown for the cut Mpp < 2m + 15 MeV. The total 
number of the pp7r~ events is reduced approximately by a factor 3: We obtain ANt ~ 2000 
triplets, in the excellent agreement with the estimate of Eq. p^). The RMS deviation 
cr ~ 65 is somewhat larger. One can expect that the a is reduced by a factor a/3 = 1.7, 
while we get only a reduction by a factor of ~ 1.3. This deviation can, however, be explained 
by statistical fluctuations due to the relatively small number of the bins (= 10). 

The angular cut | cos| < 0.75 reduces the background by a factor of two, while 
the number of the detected d' events is reduced by a factor 3/4. The results for the double 
cut Mpp < 2m + 15 MeV and | cos| < 0.75 are shown in Fig. 6 (c). Now we get 
only ANt ~ 1100 triplets. The RMS deviation a ~ 30 is again in good agreement with the 
Poisson law, i.e. cr ~ y/ANt ~ 33. 

To summarise, the selection of the pp7i~ events suppresses the background by a factor 
12. The second cut over the invariant masses of the two protons suppresses the background 
by a factor 3. The third cut over the proton scattering angles suppresses the background by 
a factor 2. The number of the detected d' events decreases by a factor 1/3 in the hrst case, 
decreases by a factor 2/3 in the second case, and decreases by a factor 3/4 in the third case. 
The number Nc, which is the ratio between the background Ant and square of the signal 
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n% (see Eq. (0). decreases by a factor 1/2. 

The fluctuations of the pp'K~ background summed up over all ion collisions are described 
reasonably well by the Poisson statistics, and so Eq. (p2|), on which we base the estimate 
for the total number of collisions needed to reveal the d' peak, are apparently justihed. 

At the same time, fluctuations of the pp'K~ background from one collision to another 
due to the variations of the impact parameter and due to variations of the coordinates and 
momenta of the nucleons inside of the colliding nuclei, are significant. These fluctuations 
are analysed in Fig. 7 with no cuts (a), with the cut Mpp < 2m + 15 MeV (b), and with 
the double cut Mpp < 2m + 15 MeV and | cos1 < 0.75 (c). The histograms (a) - (c) 
show the number of simulated collisions as a function of the number of all possible pp7r~ 
combinations per single collision in the bin at the d' dibaryon mass Md' = 2.063 GeV with 
width AMt = 1 MeV. The reduction of the background by a factor 3 on the histogram 
(b) and further by a factor 2 on the histogram (c) as compared to the histograms (a) and 
(b), respectively, is clearly seen, while the RMS deviations are 2-3 times greater than one 
could expect from the Poisson distribution. It means that the distinction between different 
collisions has essentially a dynamical origin. If all triplets from all collisions are mixed, the 
Poisson distribution starts to work, as we discussed above. 

With all the cuts taken into account, the number of the ppTT~ events in a bin AMt = 1 
MeV around the d' mass per one collision equals An* = 11 (see Fig. 7 (c)). This number 
should be compared with the number of the detected d' dibaryons in the channel ppTi~ 
Hd' ~ 0.02 ■ (1/3) ■ (2/3) ■ (3/4) = 0.0033. It means that the number of the Au + Au 
collisions needed to reveal the d' peak at a 6cr level is 

^ 3.6 10^. 30 

The A dependence of the number of the background events and of the ratio Arit/n^i are 
analysed in Figs. 3 (b) and (c). The results of the QMD simulations confirm essentially the 
naive estimates of Sect. 3: The value Ant increases as A^, while the ratio Ant/n^, increases 
linearly with A. In the general case, Eq. (30) reads 
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Nc>2 10 ^ ■ A. 


( 31 ) 


The limit on the value Nc is inverse proportional to square of the d' production cross 
section and depends linearly on the interval AM* = max(rrf/, where [AM^Y^p is 

the resolution of the detectors. The d' cross section in Eq. m is well dehned provided that 
the signal seen at CELSIUS [Q is really connected to the d' dibaryon. Otherwise the used 
value of app^d'TT+ gives the upper limit for the d' production cross section, for the number 
rid' of the detected d' dibaryons, and the minimum value for the number Nc- 

VI. SUMMARY 

We considered the production of the d' dibaryons in heavy ion collisions. Due to smallness 
of the cross section for the d' dibaryon production in the elementary NN collisions, the 
production of the d' dibaryons in heavy ion collisions was considered as a perturbation to 
the reaction dynamics. The number of the d' dibaryons was calculated using the QMD model 
for and ^^'^Au central collisions at the energy 1000 MeV per nucleon. We 

obtained rid' ~ 10“^ ■ A dibaryons per single collision. This number was compared with the 
background in the invariant mass distribution of the ppTi~ events. The comparison showed 
that at least 2 10®-A collisions are needed to make the d' signal statistically significant. Light 
nuclei are more promising candidates for the d' searches. If data on the ppTi~ invariant mass 
distribution from earlier experiments are already available, the search for the d' peak can be 
made both with light and heavy nuclei. This is a realistic task to collect statistics sufficient 
to reveal the d' peak, and thus the search for the d' dibaryon seems to be a possibility which 
should be probed first by analysing the available heavy ion data accordingly. 
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APPENDIX. THE ENHANCEMENT FACTOR rjp FOR REACTION D' NNtt IN 
THE EFFECTIVE-RANGE APPROXIMATION 


The Jost function in the effective range approximation has the form (see [17], Ch.9) 


1/D(k) = 


k + ia 
k + i'y 


where 


1 / N 1 

-r^[a - 7 ) = 1 , -r^a-f = 1/a. 


Here, a is the scattering length, fm the effective radius, and k is the nucleon momentum 
in the center-of-mass frame of two nucleons. In our case, a = 23.7 fm, re = 2.67 fm, and 


A; = y - TV). 

In Fig. 1, the value /3{T) = ^ln\D{0)/D{T)\‘^ where T = 2k‘^/m is plotted (dashed 


curve). 


The factor rjp is obtained by averaging the value \1/D{k)\^ over the phase space of the 


NNti system 


7]f = 




/ - T^)dT^ 

This expression can be transformed as follows 

TT Jc k^z + ^ 

The contour integral runs around the cut [0,1], Re(z) decreases when Im{z) is positive and 
increases when Im{z) is negative. Here, 

We move the contour C to inhnity and obtain 

riF = l + 4(a^ - - if/k^. 

Substituting a zz 2/re >> 7 ~ l/o and ko = 210 MeV >> 7 , we get 

rjF ^ 1 + lQ/{rekof ^ 3. 

For T = 800 MeV and k = JhmT = 600 MeV, 7 / = \1/D{k)\‘^ 1. 


18 



REFERENCES 


[1] R. Bilger et al. Phys. Lett. B269 (1991) 247. 

[2] R. Bilger, H.A. Clement and M.G. Schepkin, Phys. Rev. Lett. 71 (1993) 42. 

[3] M.G. Schepkin, O. Zaboronsky and H. Clement, Z. Phys. A345 (1993) 407. 

[4] L.Ya. Glozman, A. Buchmann and A. Faessler, J. Phys. G20 (1994) L49. 

[5] H. Clement, M.G. Schepkin, G.J. Wagner and O. Zaboronsky, Phys. Lett. B337 (1994) 
43. 

[6] G.J. Wagner, Prog. Part. Nucl. Phys. 34 (1994) 345 . 

[7] G. Wagner, L.Ya. Glozman, A.J. Buchmann and A. Faessler, Nucl. Phys. A594 (1995) 
263. 

[8] W. Brodowski et al. Z. Phys. A355 (1996) 5. 

[9] A. Faessler, A.J. Buchmann, M.I. Krivoruchenko and B.V. Martemyanov, Phys. Lett. 
B391 (1997) 129; J. Phys. G24 (1998) 791. 

[10] A. Faessler, A.J. Buchmann and M.I. Krivoruchenko, Phys.Rev. C56 (1997) 1576; C57 
(1998) 1458. 

[11] A.J. Buchmann, A. Faessler and M.I. Krivoruchenko, Ann. Phys. 254 (1997) 109. 

[12] J. Aichelin, Phys. Reports 202 (1991) 233. 

[13] J. Aichelin and H. Stocker, Phys. Lett. B176 (1986) 14. 

[14] C. Fuchs, L. Sehn, E. Lehmann, J. Zipprich and A. Faessler, Phys. Rev. C55 (1997) 
411. 

[15] V.S. Uma Maheswari, C. Fuchs, A. Faessler, L. Sehn, D. Kosov and Z.S. Wang, Nucl. 
Phys. A628 (9198) 669. 


19 



[16] S.M. Kiselev, Nucl.Phys. A579 (1994) 643. 

[17] M.L. Goldberger and K.M. Watson. Collision Theory, John Wiley and Sons, NY (1965). 

[18] Richard A. Arndt et al, Phys. Rev. D28 (1983) 97. 

[19] Frank Tabakin, Phys. Rev. 174 (1968) 1208. 

[20] J.W. Harris et al., Phys. Rev. Lett. 58 (1987) 463. 

[21] D. Pelte and the FOPI Colloboration, Z. Phys. A357 (1997) 215. 


20 



FIGURE CAPTIONS 


Fig. 1: The function /3(T) = ^lYL{ri{T)/ri{0)) versus the laboratory kinetic energy T of the 
incoming nucleon. The solid curve represents the calculation with the Jost function 
(7), while the dashed curve is due to the effective-range approximation described in 
Appendix. 

Fig. 2: The cross section app^d'iT+ versus the laboratory kinetic energy T from the d'-K threshold 
up to 1000 MeV for the reaction pp . 

Fig. 3: (a) The d' production rate n^/ versus the atomic number A of the colliding nuclei 

and Au at 1 A-GeV incident energy. The errors display variations of 
the value Ud' due to variations of the impact parameter and initial coordinates and 
momenta of nucleons in the colliding nuclei. The solid straight line is the linear fit 
of the value n^/. (b) The number of the pp'K~ triplets per 1 MeV of the ppn~ 
invariant mass in vicinity of the d' peak versus atomic number A of the colliding nuclei. 
The errors give statistical fluctuations. The solid curve is the A^ £t of the value Arit. 
(c) The ratio between the background ppn~ events and square of the number of the d' 
dibaryons versus atomic number A. The solid straight line is the linear fit of the ratio 
Aritlnj,. 

Fig. 4: The proton polar angular distribution in the center-of-mass system of the colliding 
nuclei in central (6 < 3 fm) ^^'^Au -|- ^^'^Au collisions at 1 A-GeV incident energy. 

Fig. 5: The distribution of the ppn~ background events in central (6 < 3 fm) 

collisions at 1 A-GeV incident energy, (a) The pp7T~ background in the interval 2-^3 
GeV of the invariant mass Mt for a single typical collision, (b) The pp7i~ background 
near threshold averaged over Nq = 100 collisions. According to Eq. ([T6|) , the distri¬ 
bution has a simple threshold behaviour: ANt = C{Mt — MthY- It is seen that the 
parabolic curve is an excellent fit to the distribution near the threshold, (c) The total 
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pp'K~ background at threshold (see scale of Mt) for Nc = 100 collisions. The straight 
solid line gives the best £t for these 10 bins. The deviations of the AiV^ from the mean 
values indicated by the solid line are in good agreement with the Poisson distribution 
(see Figs. 6 and discussion in the text). The value AiVt = {dNt/dMt)AMt is the 
number of the pp7r~ events in the interval AMt. In the histograms (a), (b), and (c), 
AMt = 10, 1, and 1 MeV, respectively. 

Fig. 6 : The nnmbers of the bins of Fig. 5 (c) with the fixed deviations from the mean valnes 
of the ANt, indicated by the solid line, with no cnts (a), with the cnt over the two- 
proton invariant masses Mpp < 2m + 15 MeV (b), and with the double cut over the 
two-proton invariant masses Mpp < 2m -|- 15 MeV and the proton scattering angles 
I cos0p™'^'| < 0.75 in the center-of-mass system of the colliding nuclei (c). The RMS 
deviations a = 80, 65, and 30 should be compared with the mean values ANt = 5500, 
2000 , and 1100 in the cases (a), (b), and (c), respectively (the mean values ANt/Nc 
with Nc = 100 can be read off from the upper right corners of Figs. 7 (a) - (c)). The 
flnctuations of the backgronnd a ~ ANt do not contradict to the Poisson statistics. 

Fig. 7: The distribution of 100 simulations of Au + Au collisions as a function of the 
nnmber of the ppTT~ combinations (triplets) per single collision in the region Md' ± 0.5 
MeV with no cnts (a), with the cut over the two-proton invariant mass Mpp < 2m+ 15 
MeV (b), and with the donble cut over the two-proton invariant mass Mpp < 2m + 15 
MeV and the proton scattering angles | cos0p™'^'| < 0.75 in the center-of- mass system 
of the colliding nuclei (c). 
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